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We investigate the QCD magnetic susceptibility Xq ^or flavor SU(2) at finite temperature (T) 
beyond the chiral limit, using the liquid instanton model, defined in Euclidean space and modified 
by the T-dependent caloron solution. The background electromagnetic fields are induced to the 
QCD vacuum, employing the Schwinger method. We first compute the scalar (chiral) and tensor 
condensates as functions of T as well the current-quark mass m, signaling the correct universal chiral 
restoration patterns. It turns out that Xq^ given by the ratio of the two condensates, is a smoothly 
decreasing function of T, showing about 20% reduction of its strength at the chiral transition T = To, 
in comparison to that at T = 0, and decreases almost linearly beyond To for m j^ 0. We observe 
that the present numerical results are in qualitatively good agreement with other theoretical results, 
including the lattice simulations. Finally, we examine the effects of the external magnetic field on 

C^ the tensor-polarization VEV, resulting in that it plays the role of the chiral order parameter. 
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^ I. INTRODUCTION 

f^ The low-energy quantum chromodynamics (QCD) manifesting nonperturbative natures have been investigated 

Qh extensively in many different ways, such as the lattice QCD (LQCD), effective QCD-like models, QCD sum rule 

f-^ (QCDSR), and so on. To understand the nonperturbative features of QCD, it is necessary to scrutinize the QCD 

q^ vacuum structure, which governs the breakdown of relevant symmetries of QCD and the phase transitions of the QCD 

,^ matter. Note that, recently, the electromagnetic (EM) properties of the QCD vacuum at finite temperature (T) and/or 

'~~' finite quark chemical potential (/z^) have attracted much attention from experiments [1] as well as theories [IHl], being 

I together with the energetic progresses of the heavy-ion collision experiments. Among the relevant physical quantities, 

J>. which are sensitive to the U(l) EM interactions, the QCD magnetic susceptibility is one of the important ones, due 

l/~j to its impact in theories and experiments. The QCD magnetic susceptibility for a quark flavor q, Xq is stands for a 

^>0 response of the scalar (chiral) condensate to the external EM background field, and defined in terms of the vacuum 

CN expectation value (VEV) for the tensor-polarization operator (TP-VEV), {q^ cr^^^q)^^, in Euclidean space iSnBj: 

"^ {q^cf^i,q)-E-M = e.qF^^{iq^q)Xq, (1) 

o 

C^ where e^ and F^^ denote the quark electric charge and the EM field strength tensor. The subscript EM in the 
^^ left-hand-side of Eq. (fTl) stands for the existence of the external EM field. Note that TP-VEV is a linear function of 
^ F^^ in the leading order. The magnetic susceptibility is also pertinent to the photon distribution amplitude |7| as 
well as the anisotropy in terms of the para- or dia-magnetism [8] . This physical quantity was studied up to now in 
QCDSR 911, effective quark models [3 [12], holographic QCD (hQCD) 13, 14], operator product expansion with 
S the pion dominance (OPE-I-PD) [T5], and LQCD [TMT5] . and so on. Note that, in the previous work, we computed 
Xq using the liquid instanton model (LIM) [121 120] at finite density [H], in which we noticed that x^ > with 
the Euclidean metric, indicating the diamagnetic matter [xq < for the Minkowski metric). Taking the present 
heavy-ion collision experiments at low quark density into account, in the present work, we are focusing on the QCD 
magnetic susceptibility for the flavor SU(2), SU(2/), at finite T and fiq = 0, using LIM, modified by the T-dependent 
caloron solution [321 HI]- We mention that all the calculations are performed in Euclidean space. The relevant 
instanton parameters, such as the average (anti)instanton size p and inter-(anti)instanton distance R are modified as 
functions of temperature, employing the trivial-holonomy (Harrington-Shepard) caloron solution [ 2Hl26| . As a result, 
these modifications show the partial restorations of the spontaneous breakdown of the chiral symmetry (SBxS) as 
T increases. Relating to those model parameters, it is worth mentioning that our model (renormalization) scale is 
taken as /i ~ \/2/p « 0.85 GeV at T = [22|. The T-dependent effective quark mass is computed numerically as 
an order parameter for SBxS, employing the T-modified effective thermodynamic potential [26) . By doing that, we 
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observe correct universal chiral restoration patterns for the zero and finite current-quark masses, i.e. the second-order 
and crossover chiral restoration patterns. Being equipped with those ingredients, we present the numerical results 
for the chiral condensates (iq^q), tensor condensate (q^Eg), their ratio corresponding to the magnetic susceptibility 
{q'^T.q) / {iq^ q) = Xq^ ^^d TP-VEV {q^ (J^vq)-EM- Here, we write the definition of the tensor condensate as follows: 

{q'^q) ^ i^!^^^ = x,{^<i\) for F,. ^ 0. (2) 

For instance, in Ref. [IH], the authors defined the tensor condensate as (q^'Sq) = — t^. 

From the numerical calculations, it turns out that the chiral condensate exhibits the correct chiral restoration 
patterns as expected from those of the effective quark mass as mentioned above. The chiral transition T is given 
as To — (166, 170) MeV for the zero and finite current-quark masses, i.e. m — 0, and 5 MeV, considering the 
SU(2j) flavor symmetry. Note that the values for the chiral condensates are obtained as (iq^q) = (258, 260MeV)^ 
for m — (0, 5) MeV at T = 0. These values are well compatible with its empirical one (240 ± lOMeV)'^. Using the 
same parameter sets, we also compute the tensor condensates, which also manifest the universal chiral restoration 
patterns. At zero temperature, we have {q^'Sq) — (52,49) MeV for m — (0,5) MeV at T = 0. Interestingly, we 
observe bump structures in the tensor condensate curves for m ^ a.t T = (50 ^ 60) MeV, due to the nontrivial 
interference between the constituent- and current-quark masses. These values are again well-compatible with the 
known theoretical estimations O [7l [T0Ul2l [T2j [T8] . When we compare the temperature-dependent behaviors of the 
present numerical results with the lattice data [TH] , extrapolated to the physical u-quark mass, at the renormalization 
scale fJ, = 1 GeV, we observe qualitatively good agreement with them, but sizable deviations also appear in the vicinity 
of the chiral transition temperature Tq. At the chiral phase transition T, Tq = 170 MeV, the tensor condensate for 
TO 7^ becomes 21 MeV, whereas it is zero for the chiral limit. 

It turns out that the magnetic susceptibility decrease steadily with respect to T. We obtain their typical values, 
estimated as Xq = (3.03, 2.77) GeV~^ for to, = (0, 5) MeV at T = 0. Again, these values are well matched with those 
from the LQCD and other effective models. Beyond Tq = 170 MeV, the curve for the magnetic susceptibility behave 
almost as a linearly decreasing one as T increases. At Tq, we observe about 20% decreases in their strengths, in 
comparison to those at T = 0. TP-VEV is a linear function of the external magnetic field (e^S) in the leading order. 
The slope of the TP-VEV line with respect to CqB decreases as T increases, since the tensor condensate plays the 
role for its slope value, signaling the (partial) restoration of SB^S. Thus, TP-VEV can be considered as a chiral order 
parameter. Consequently, we observe that it vanishes at Tq = 166 MeV in the chiral limit, whereas remains finite 
beyond Tq for to ^ 0, because of the crossover chiral phase transition. 

We organize the present work as follows: In Section II, we bricfiy introduce the liquid instanton model (LIM) and 
how to compute the magnetic susceptibility in terms of the field theoretical manner. In Section III, the temperature 
modifications of the relevant model parameters are performed using the trivial caloron solution. We also show the 
correct universal chiral restoration patterns, computed within the present model. The numerical results for the chiral 
and tensor condensates as functions of temperature are presented with relevant discussions in Section IV. In addition, 
the magnetic susceptibility is estimated and compared with other theoretical estimations. Final Section is devoted to 
summary, conclusion, and future perspectives. 

II. EFFECTIVE ACTION VIA THE INSTANTON- VACUUM CONFIGURATION 

In this Section, we introduce the liquid-instanton model (LIM) briefly as a theoretical framework to study the 
magnetic susceptibility. Details on the present framework can be found in Refs. (THl US]. The effective action for 
SU(2/) via the instanton vacuum can be written in Euclidean momentum as follows [5l 119): 



5eff [to, A^,T^,] = -Sp,j^^ In [ip + ^TO + iM{d^) + a ■ T] , (3) 

where Sp^ t represent the functional trace running over the color (c), flavor (/), and Lorentz index (7). The U(l) 
covariant derivative reads iD^ = idf^ + EqA^, in which Cq stands for the electric charge of a quark, rh stands for the 
current-quark mass matrix, d\a,g{mu,md). Throughout this work, we assume the SU(2^) symmetry for the quark 
masses, i.e. to„ ss rud « to. = 5 MeV for the cases beyond the chiral limit. The effective quark mass M(9^) is 
generated from the nontrivial interactions between the quarks and (anti)instanton via the quark zero mode |19j . and 
it reads 



M(a2) = MQF^{d^) = Mq 



2 

2 + p2|a|2 



2 



M{k^) = Mq 



2 

2+p2fc2 



2 



= Mfe. (4) 



Here, p denotes the average (anti)instanton size ~ 1/3 fm [TO], and Mq indicates the constituent-quark mass at zero 
virtuahty. In the last step of Eq. m , we wrote the effective mass in the Euchdean momentum space. We defined the 
antisymmetric tensor a^y = i(7p7i/ — 7i/7p)/2 with the external tensor source field T^"^. From the effective action in 
Eq. (pl, one can write the chiral condensate space by performing the functional derivative of the effective action with 
respect to m as follows: 



{iq^q) = AN, 



m 



m 



fc2 



(5) 



where we assign as J^, = / -A^ for convenience. We also have used a simplified notation M^ — ra-\- Mk- Employing 
the phenomenological values for the model parameters 1/p ~ 600 MeV and Mq w 350 MeV, which are chosen to 
reproduce the pion weak-decay constant F^^ « 93 MeV, we obtain (iq^q) ~ (250 MeV)"^ in the chiral limit. Note that 
this value is well compatible with its empirical values j27' . 

Similarly, the matrix element in the left-hand-side of Eq. (fTl), i.e. TP-VEV can be evaluated by performing the 
functional derivative with respect to T^'^ in the presence of the EM background field, induced by the Schwinger 
method [SmilllH]: 



{q^'^p^vq)i 



1 95eff[m,A^,rn 



Nf 



STA"^ 



T=0 



(6) 



By expanding Eq. (6| in terms of e^, we obtain the following expression in the leading order ex ©(cg), according to 
e^ < 1 for n > 2 U^- 



{q^(^^iyq)^M = '^Nc{eqF^^) 



(2^ 



Mk + Nk m 

(fc2 + 1/2)2 (fc2 + ^2)2 



e„F, 



gJ- ti,v^ 



?^Sg). 



(7) 



Here, we assign {q^I]q) as a tensor condensate for convenience as mentioned previously. Here, E stands for a scalar 
operator to satisfy Eq. ([7| , and its analytic form does not make any impact on the final results of the present work. 
It is worth mentioning that the tensor condensate is defined alternatively as (q^'Sq) = —Tq with the Minkowski 
metric [18]. The mass-derivative term Nk is defined as 



Nk = -k 



.dMk SMoip^k' 



2l,2\ 



afc2 (2 + p2/j2)3 



(8) 



Notice that TP-VEV in Eq. (I?]) is a linear function of the field-strength tensor Ffj_i, in the leading order expansion of 
Cq. This observation is consistent with other theoretical studies '^. Details of the derivation of Eqs. JtI) and (IS]) can 
be found in our previous work (21j and references therein. Substituting Eqs. ([7]) and ([5]) into Eq. (JlV we have the 
following equation for the QCD magnetic susceptibility as follows: 



Xq 



{iq^q) 



Mk + Nk 
(fc2 + Af 2)2 



(/c2 



Mk 



fc2 



m 



(9) 



III. TEMPERATURE-DEPENDENT EFFECTIVE QUARK MASS FOR FLAVOR SU(2/) 

In this Section, we would like to briefly discuss how to compute the effective quark mass Mq in Eq. Q as a function 
of T, and to develop the T-dependences for the model parameters. In Refs. [25] j we derived it by using the caloron 
distribution with trivial holonomy, i.e. Harrington-Shepard caloron |24[ 125] . Firstly, we want to explain briefly how 
to modify p and R as functions of T, using the caloron solution. Details can be found in Ref. [26 . An instanton 
distribution function for arbitrary N, and Nf can be written with a Gaussian suppression factor as a function of T 
and an arbitrary instanton size p for pure-glue QCD |25j : 



d{p,T) = C_ 



N, 






/3^^ p'- 



•exp[-{ANT^+Pjnp^)p^] 



(10) 



We note that the CP-invariant vacuum was taken into account in Eq. (10 1, and we assumed the same analytical 



form of the distribution function for both the instanton and anti-instanton. Note that the instanton number density 



(packing fraction) N /V 



1/i?^ and p have been taken into account as functions of T implicitly. For simplicity, 



we take the numbers of the anti-instanton and instanton are the same, i.e. Nj = Nj = N . We also assigned the 
constant factor in the right-hand-side of the above equation as C for simplicity. The abbreviated notations are also 
given as: 



An, = 



51n[ARsPcut], /3 
11 



6 



-Nr-l 



-61n[AR,s(i?)], Cn, 

27 r N, 
4 m -I 



4.60 g-'^-^^O'^sNc 

7r2(7Ve-2)!(iV,-l)!' 



6 = 



3 ' 



(11) 



Note that we defined the one-loop inverse charge /3 and /? at certain phenomenological cutoff pcut and (R) « R. 
Ars denotes a scale, depending on a r enor malization scheme, whereas V3 for the three-dimensional volume. Using the 
instanton distribution function in Eq. ( 10 1, we can compute the average value of the instanton size p^ straightforwardly 
as follows [25;: 



p'iT) 
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fdpd{p,T) 
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where v ~ [h — 4)/2. It can be easily shown that Eq. (12) satisfies the following asymptotic behaviors 
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Here, the second relation of E q. ([l3| indicates a correct scale-temperature behavior at high T, i.e., 1/p « A c>c T. 
Substituting Eq. ( 12 ) into Eq. (|10[), the caloron distribution function can be evaluated further: 



1 



j-(r) = ^AnT' 



1 



-A'nT 



ujS^n 



d(p,r)=Cp^-5exp[-J-(T)p2] 
The instanton packing fraction n can be computed self-consistently, using the following equation: 



(14) 



(15) 



where we replaced NT/V3 — » n, and r(t^ ) stands for the F-function with an argument v. Note that C and /? can be 
determined easily using Eqs. ( 12 ) and (Il5|, incorporating the vacuum values for n « (200 MeV)'* and p « (600 MeV)~^: 
C w 9.81 X 10"'* and /3 « 9.19. Finally, in order for estimating the T-dependence of Mq, one needs to consider the 
normalized distribution function, defined as follows: 



dN{p,T) 



Jdpd{p,T) 



>-5 



■F-(T)exp[-.F(T)p2] 
F(z.) 



(16) 



Here, the subscript N denotes the normalized distribution. For brevity, we want to employ the large-TVc limit to 
simplify the expression for dM{p,T). In this limit, as understood from Eq. (16 1, d^i^p^T) can be approximated as a 
(5-function: 



N, 



lim dNip,T)^d[p~p{T)]. 



(17) 



The numerical result for the trajectories of p{T) (solid) and 1/R{T) (dot) are given in the panel (a) in Figure ^ 
Here we choose p{0) ~ 1/3 fm and R^ 1 im for all the numerical calculations. These values are phenomenologically 
preferred in the present model '19]. The curve for p{T) shows that the average (anti)instanton size smoothly decreases 
with respect to T, indicating that the instanton ensemble gets diluted and the nonperturbative effects via the quark- 
instanton interactions are diminished. At T = (150 ~ 200) MeV, which is close to the chiral phase transition T, the 
instanton size decreases by about (10 ~ 20)% in comparison to its vacuum value. Considering that the instanton 
size corresponds to the scale parameter of the model, i.e. UV cutoff mass, p « l/{\/2p), the T-dcpendcnt cutoff 
mass is a clearly distinctive feature in comparison to other low-energy effective models, such as the NJL model. In 
addition, we also show the T dependence of the average (anti)instanton number density or (anti)instanton packing 
fraction, N/V ~ ^/R'^, in the panel (a) of Figure^ Again, the instanton number density get decreased as temperature 
increases: The instanton ensemble gets diluted. We will use these two temperature-dependent quantities for computing 
the chiral and tensor condensates, and TP-VEV in the next Section. 
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FIG. 1: Average (anti)mstanton size p ~ 1/A [fm] and (anti)instanton packing fraction {N/VY''^ [GeV] as functions of T, 
computed from the Har ring ton- Shepard caloron distribution 24, 25j in the panel (a). Effective quark mass at zero virtuality, 
Mo computed from Eq. ( 19 1 as functions of T for ra = Q (solid) and m = 5 MeV (dot), signaling the second-order and crossover 
chiral phase transitions, respectively, in the panel (b). The vertical lines indicate the chiral-phase-transition temperatures 
To = (166, 170) MeV for m = (0, 5) MeV. 



Now, we are in a position to discuss the T dependence of the constituent quark mass Mq as a chiral order parameter. 
As in Ref. [5S], the LIM thermodynamic potential per volume in the leading large- A^c contributions at zero quark 
chemical potential can be written as follows: 



n 



LIM 



A^ 
V 



1 - In ■ 



A^ 



20-2 - 2NcNf 



(Pk 

(2^ 



2Tln 



1 



(18) 



where A and M represent a Lagrange multiplier to exponentiate the effective quark-instanton action and an arbitrary 
massive parameter to make the argument for the logarithm dimensionless. a stands for the VEV for the isosinglet 
scalar meson field corresponding to the effective quark mass. The quark energy is defined by E^ — k'^ + M|. In the 
leading large- A^c contributions, we have the relation 2a^ = N/V |26|- The gap equation can be derived from Eq. (18) 
by differentiating JIlim by the Lagrange multiplier A: 
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Note that one can write the instanton packing fraction in terms of the effective quark mass Mq and p |19) : 



A^ 
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CoNeMl 



(19) 



(20) 



The value of Cq locates in (1/3 - 1/4) for 1/p « 600 MeV, Mq « (300 ~ 400) MeV and N/V w (200 - 260 MeV)^ for 
vacuum [3D]. We choose Co = 0.27 to reproduce Mq = (340 ~ 350) MeV at (T, /x) = in the chiral limit. After solving 
Eq. ( [l9| ) with respect to Mq numerically, the numerical results for Mq as a function of T are given in the panel (b) 
of Figure [I] for the zero and finite current quark mass: m — (solid) and to = 5 MeV (dot). These results indicate 
correct universal patterns for the chiral phase transition like the those of the Ising model, i.e. the second-order chiral 
phase transition for the massless fermion and the crossover for the finite mass. From those numerical results, the 
phase transition T for the two chiral restorations are obtained as Tq sa (166, 170) MeV for m = (0, 5) MeV. The 
transition temperatures are indicated by the thin solid vertical lines in the panel (b) of Figure [I] Detailed discussions 
for the chiral phase structure within the present model is given in Ref. 31j. 

To evaluate Xq as a function of T, we redefine Eq. ^ with the fermionic Matsubara formula: The integral over the 
fourth momentum, k^, is compactified into a summation over the Matsubara frequency: 



Xg 
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(21) 



where the fermionic Matsubara frequency Wn = (2n + l)7rT and the three-dimensional integral is given in a simplified 



Here, we use the notation En = k^ 



whereas k denotes the three momentum of the 



notation J, = j /„ ^3 . HCIC, WC use UHC UULClLiUii JJyg 

quark. Here is one caveat: Introducing the Matsubara frequency, we assume that the effective quark mass in Eq. (HI) 
is simplified by fc4 — > 0: 



Mk -^ Mk = Mo 



2 + p2fe2 



(22) 



and the same for Nk — > N^, in Eq. ([8]). We have verified that this simplification makes the problem in hand much 
convenient and simplified for the analytic as well as the numerical calculations, and does not make significant deviations 
from the full calculations, as shown in many successful applications |221 1231 1^ . Then, the summations over Wn in 
Eq. (21) can be analytically performed, and we defined the following functions: 



f2(E) = ^"--scch^ 



2T 



Tsinh 



-E 



1 f E 

ME) = — tanh — 



Using Eq. (23), Eq. (21 ) can be then rewritten finally as 



Xq 



{iq^q) 



j [{Mk + Nk)f2{Ek) - mh{E„)] I jy [Affc/i(^fc) - m/i(So)] 



(23) 



(24) 



Note that all the quantities, Xq, and chiral and tensor condensates, are positive-real valued functions of T and m in 
the present work with the Euclidean metric. Hence, our theoretical results exhibit the diamagnetism, considering that 
X^ < for the Minkowski (M) metric. 



IV. NUMERICAL RESULTS AND DISCUSSIONS 



In this Section, we demonstrate the numerical results with relevant discussions. First, we show them for the 
chiral condensates (iq^q) as functions of T in the panel (a) of Figure [2] The (solid, dash) curves correspond to 
the condensates with m = (0, 5) MeV. In what follows, numerical values given in the form of {x, y) represents the 
theoretical results for m, = (0, 5) MeV, respectively, unless otherwise stated. The vertical straight lines denote the 
chiral transition T, Tq = (166, 170) MeV. At T = 0, we observe (iq'^q)^/^ = (258 MeV, 260 MeV) at T = 0. Here, we 
only observe small deviations depending on m. These numerical values are slightly larger than the empirical values 
(iq^q) = (240 ± lOMeV)'^, but still in qualitative agreement with them. As for the curve for m = 0, it shows the 
second-order chiral phase transition as expected and understood by Eq. ([5]), which is proportional to Mq. On the 
contrary, m becomes finite, the curves manifest the crossover transition, satisfying the universal class pattern of the 
chiral restoration. At Tq — 170 MeV, we have {iq^q) ~ (209 MeV)^ for m = 5 MeV, showing about 20% reduction, in 
comparison to that for T = 0. 

In the panel (b) of Figure pi we depict the tensor condensate (g^Sg) = —Tq with the same manner of the panel (a). 
Again, the condensates show the proper chiral restoration patterns depending on m. Interestingly, there appear bump 
structures in the curves for m ^ at T = (50 ~ 60) MeV, due to the nontrivial interference between the constituent- 
and current-quark mass terms in Eq. (ITJ). We find (q^'Sq) — (52,49) MeV at T = 0. From these values, we conclude 
that the tensor condensate decreases with respect to m which is consistent with the observation of Ref. [32]. Here, 
we want to mention other theoretical estimations for the tensor condensate for SU(2j) for T — 0. Using QCDSR 
techniques, it was studied in Refs. [TIUHIIII], which estimated it as (40 ~ 70) MeV, depending on the different 
renormalization scales fi = 0.5 GeV or /x = 1 GeV. In Ref. 6^, employing the same instanton model for vacuum, with 
slightly different model parameters, the authors calculated it in the chiral limit, resulting in (45 ~ 50) MeV, which 
is well compatible with ours by construction. Employing the NJL model and quark model (QM), it was estimated 
as 69MeV and 65 MeV, respectively, in Ref. [12| which is about 10% larger than ours. From the quenched LQCD 
simulations for SU(2c) [H] and SU(3c) [H], they observed 46 MeV and - 52 MeV at ^ = 2 GeV. A fuU SU(3c) lattice 
simulation was performed in Ref. [18], and provided (38.9 ^ 40.7) for the chiral limit and physical current-quark 
mass at ij, — 2 GeV. Considering that the renormalization constant for the tensor condensate with the running scale 
yU = (1 ~ 2) GeV is close to unity as will be shown below, our present estimations for the tensor condensates are well 
compatible with other theoretical results. The comparisons with other studies are summarized in Table |l| 

Now, we are in a position to discuss the T dependence of the tensor condensate. Bali et at estimated the T 
dependence of the tensor condensate via the full SU(3c) lattice simulation with tree-level Symanzik improved gauge 
action [18_ . They also performed renormalization-group analyses, resulting in the renormalization constants for 
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FIG. 2: (Color online) (a) Chiral condensate {iq^q) = —{qq) in Eq. ( 24 1 for m — (0, 5) MeV in the (solid, dash) lines, respectively. 
The vertical lines indicate the chiral phase transition T, To = (166, 170) MeV for {m — 0,m j^ 0). The horizontal shaded area 
denotes the range of {iq^q) — (250 ~ 260 MeV)'', which corresponds to its empirical value, (b) Tensor condensate (g^Eg) = — r. 



in Eq. ( |24[ ), represented in the same manner with the panel (a). The lattice QCD data at the renormalization scale fJ. = 1 GeV 
are taken from Ref. [18], and indicates the continuum extrapolation for the u quark with the errors containing all statistical 
and systematic errors. The data for the d quark are still within the u-quark errors. 



the chiral (scalar) (S) and tensor (T) condensates for the running scale /i 

Xq = {q^^q)/{iqU)- 
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1) GeV as follows, considering 



(25) 



Note that in their lattice simulation, they renormalized all the quantities at /i = 2 GeV with (m^ + rn^) « 6.94 MeV, 
which gives (qq) = — (269 MeV)"^. Using Z|j^ in Eq. (25), this chiral condensate value becomes — (245MeV)'^ at /i = 1 

GeV. Since our renormalization scale is about 0.85 GeV ~ V^/p, for appropriate comparison, one needs to evolve their 
data at /x = 2 GeV to those at ^ = 1 GeV. After the scale evolution, just multiplying Z^ to the data, the continuum 
extrapolated data for the u quark are given with the open dot with the sum of the statistical and systematic errors 
in the panel (b) of Figure [2] Note that the lattice data locates between the curves for m = and m = 5 MeV. 
Although there appear quantitative differences between the lattice data and ours, the overall tendency and strength 
are qualitatively comparable. Note that the chiral transition T was given Tq sa 162 MeV in the lattice simulation 
which is a only few percent smaller than ours Tq = 170 MeV. At the chiral phase transition T, Tq = 170 MeV, the 
tensor condensate for m — 5 MeV becomes 21 MeV from our calculations, showing about two-times reduction. 

In the panel (a) of Figure l3J we show the numerical results for Xq as functions of T for m = (0, 5) MeV in the (solid, 
dash) lines. Since Xq is the ratio of the chiral and tensor condensates as shown in Eq. (|9|, we can not define it in the 
chiral limit beyond Tq, while Xq for m j^ has finite values for T > Tq as shown in the panel (a). The typical values 
for Xq at T = are given by Xq = (3-03, 2.77) GeV"^ Note that we have Xq = (2-85 - 5.7) GeV"^ from the QCDSR 
methods [TKinilll], and these values are well compatible with our estimations. The LQCD simulations also estimated 
comparable values with ours [TMIS] as shown in Table ll] However, the hQCD [Tl] and OPE+PD [15] calculations 
showed considerably larger values for them: Xq — H-^ and 8.91, respectively. The effective quark models, such as 
NJL and QM, evaluated Xq ~ 4.3 and 5.25, depending on each models. From these observations, we can conclude 
that our model estimations are compatible qualitatively with other theoretical ones. At T = 0, we depict some LQCD 
results in the panel (a) of Figure [3J The solid square and circle are the estimations at // = 1 GeV from the full SU(3c) 
LQCD simulation for m — 5 after a proper scale evolution by multiplying Z^/Z=j=^ — 1.49. It turns out that the 
LQCD data match with the present numerical curve for m — approximately, ancl larger than that for m — 5 MeV 
by about 10%. The SU(2c) quenched LQCD data at /x = 2 GeV are also shown with the solid triangle and diamond 
for different T [16]. Although their estimation at 0.82 Tc with Tc = 313 MeV (triangle) is comparable to ours, one 
must be careful about that the simulations were done for N^, = 2 aX a. relatively larger renormalization scale without 
dynamic quarks. Note that their value at T = (diamond) is much smaller than ours as well as the SU(3c) LQCD 
simulation |18j . It was suggested that the magnetic susceptibility can be parameterized in terms of OPE+PD, using 
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FIG. 3: (Color online) (a) Magnetic susceptibility Xq ^^ functions of T for m = (0, 5) MeV in the (solid, dash) lines, respectively. 
The SU(3c) lattice QCD data at the renormalization scale [i = \ GeV are taken from Ref. [18j for the u-quark and d-quark, 
given in the solid square and circle. The SU(2c) quenched LQCD data at /x = 2 GeV are also shown with the solid triangle 
and diamond [TB]. The vertical lines indicate the chiral phase transition T, To = (166,170) MeV for m = (0,5) MeV. (b) 
Tensor-polarization VEV (TP-VEV), {q^ (JxyCi)'Eyi as functions of eqB for different T and ra values. See the text for details. 



the Gell-Mann-Oakes-Renner (GMOR) relation, as follows |15) : 



Xq 






2. 22 c^^ GeV" 



foriV, 



(26) 



where we have chosen the normalization /^ = V^^-n- with Fj^ = 92.4 MeV and Cy, stands for a positive real constant. 
We list the (average) values for c^^ for the various theory calculations in Table III Approximately, its value amounts 
to (1.0 ^ 2.0) for various calculations, whereas the OPE+PD and hQCD still give larger values than others as the 
magnetic susceptibility. 

Finally, we want to examine the external magnetic field dependence for TP-VEV in Eq. (fTl). As already mentioned, 
in the leading order of e^, TP-VEV is a linear function of the field strength tensor F^^. Choosing a certain configuration 
for Ff^^, we can write the following equation from Eq. ([l|): 



(gVahg)EM = eqB(q^T,q) 



(27) 



where the Lorentz indices a and b are understood to pick up the magnetic field from the field strength tensor in 
Euclidean space. B denotes the strength of the external magnetic field B = \B\. Considering the positive quark 
charge, u quark for instance, one can parameterize e^B as a real positive variable in the unit of GeV^. Note that 
the tensor condensate in the right-hand-side of Eq. (27) has been already computed as above. In the panel (b) of 
Figure [31 we show TP-VEV as functions of CgB for different T. The thick and thin lines denote the those for m = 
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TABLE I: Various theoretical estimations for the tensor condensate (q'Y:q) [MeV] and magnetic susceptibility Xg [GeV ^] for 
certain renormalization scale /i [GeV] and the current-quark mass [MeV] at T = 0. All the listed values are converted to those 
for Euclidean space. The notation LQCD;^ |% indicates the (full, quenched) LQCD simulations for Nc = (2, 3). c^, denotes the 
average value over possible Xg- 
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FIG. 4: (Color online) Tensor-polarization VEV (TP-VEV), {q^aabq)EM [GeV^] 
(left) and m = 5 MeV (right). See the text for details. 



and m — 5 MeV, respectively. It turns out that, as T increases, the slope of the lines decreases for the both current- 
quark masses. This behavior signals the (partial) chiral restoration, due to the decreasing of the tensor condensate, 
which plays the role of the slope of the line here. Note that this observation is consistent with that from the SU(3c) 
LQCD simulation [TB] , although the strength of the lines are different by about two times and the LQCD data show 
nonlinearity as T increases which is not shown for the present leading-order calculations. The strength difference can 
be understood by the different renormalization scales, i.e. n « 0.85 for ours and ^ — 2 GeV for the LQCD simulation. 
If we go beyond the leading order 0(eg), it is sure that there appears nonlinearity in TP-VEV. However, we will not 
discuss this issue in the present work, and leave it for the future works. Comparing the lines for the chiral limit and 
finite quark mass, the line slope get diminished much for m = with respect to T, since TP-VEV also plays the role 
of the chiral order parameter. In other words, TP-VEV in the chiral limit becomes zero for even finite CqB at the 
chiral transition T. To see this situation clearly, in Figure |4J we show TP-VEV as functions of eB as well as T for 
m — (left) and m — 5 MeV (right). As shown in the left panel, TP-VEV behaves linearly with respect to eB, and 
decreases its strength as T increases. Then, it vanishes at Tq = 166 MeV. In contrast, due to the crossover chiral 
phase transition, TP-VEV for m = 5 MeV remains finite even beyond Tq = 170 MeV. 



SUMMARY AND CONCLUSION 



In the present work, we have investigated the QCD magnetic susceptibility for SU(2/) as a function of temperature, 
beyond the chiral limit. For this purpose, we employed the liquid-instanton model (LIM), being modified by the 
trivial-holonomy caloron solution. We calculated the chiral {iq^q) and magnetic (q^^q) condensates, QCD magnetic 
susceptibility x^, and tensor-polarization VEV (TP-VEV) {q'^(Tfj_i,T,q)EM, numerically in Euclidean space. We com- 
pared our results with various theoretical estimations for those nonperturbative quantities. Important observations 
in the present work are given as follows: 



• The chiral and tensor condensates manifest the role of the chiral order parameter, showing correct chiral restora- 
tion patters. We observe (g^Sg) = (52,49) for m — (0,5) MeV at T = 0. These values are well comparable 
with the widely accepted one ~ 50 MeV. The T dependence of the tensor condensates are also compared with 
the LQCD data, and show qualitative agreement. At the chiral phase transition T, To = 170 MeV, the tensor 
condensate for m ^ Q becomes 21 MeV, whereas it is zero for the chiral limit. 

• We find that the magnetic susceptibility is a smoothly decreasing function of T. The typical values for them 
are estimated as Xq = (3.03, 2.77) GeV^^ for m = (0, 5) MeV at T = 0. Again, these values are well compatible 
with the results from the LQCD and other effective models. Beyond Tg = 170 MeV, the curve for the magnetic 
susceptibility behave almost as a linearly decreasing one with respect to T. At Tq, we observe about 20% 
decreases in their strengths, in comparison to those at T = 0. 

• TP-VEV is a linearly increasing function of the external magnetic field {cqB) in the leading order. The slope 
of the TP-VEV line with respect to tqB decreases as T increases, since the tensor condensate plays the role for 
its slope, signaling the (partial) restoration of SB^S. In that way, TP-VEV can be considered as a chiral order 
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parameter. Hence, we find that it vanislies at Tq = 166 MeV in the chiral limit, whereas it remains finite beyond 
To for 771 7^ 0, due to the crossover chiral phase transition. 

From the above observations, we can conclude that the present model calculations have revealed meaningful and 
reliable results. Since the effects of the EM fields to QCD vacuum has been one of the most energetically progressing 
objects nowadays, it is meaningful to study more various nonperturbative quantities, being sensitive to the gauge 
field, such as the mixed quark-gluon condensate {qa ■ Gq)EM, where G^i, represents the gluon field strength tensor, 
for instance. Moreover, as discussed in the previous Section, calculations beyond the linearity on the magnetic field, 
i.e. beyond the leading order, can provide interesting modifications to the present results. Related works are under 
progress and appear elsewhere. 
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